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Abstract 

We present results on the X(3872), produced in B+ X{3872)K+ and 5° X{3872)K^ decays 
where X(3872) J/V'tt+tt". We report the first statistically significant observation of B*^ 
X(3872)Kg and measure the ratio of branching fractions to be g(^^+nx(3872)i^+) = 0.82±0.22±0.05, 
consistent with unity. The mass difference between the X(3872) states produced in B^ and B^ 
decay is found to be SM = Mxk+ - Mxro = (+0.18 ± 0.89 ± 0.26) MeV/c^, consistent with zero. 
In addition, we search for the X(3872) in the decay B^ X(3872)ir+7r-, X(3872) J/^/^tt+tt-. 
We measure B{B^ X(3872)(K+^-)a//,>) x ^(X(3872) ^ J/07r+^-) = (8.1±2.0t{i) x IQ-^ and 
we set the 90% C.L. limit, B{B^ X(3872)K*(892)0) x ^(X(3872) ^ J/V'tt+tt-) < 3.4 x IQ-^. 
The analysis is based on a 605 fb~^ data sample collected at the T(4S') with the Belle detector at 
the KEKB collider. 
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INTRODUCTION 



The X(3872) was first observed in the charged i?-meson decay —>■ X{3872)K~^, 
X(3872) — > J/ip-n^-n^ by the Belle Collaboration [2]. Its existence has been confirmed by 
the CDF and DO Collaborations [2| through its inclusive production in proton-antiproton 
collisions. The discovery mode — > X(3872)i^'"'" has also been confirmed by the BaBar 
Collaboration g. The X(3872) mass, combining all measurements in this final state, is [3] 

= (3871.2 ±0.5)MeV (1) 

which is at the threshold for the production of the charmed meson pair D^D*^ . Recent 
studies from Belle and CDF that combine angular information, and kinematic properties of 
the TT^TT^ pair, strongly favor a J^*-^ = or 2 ^ assignment 0,0,8. The X(3872) does 
not appear to be a simple quark model qq meson state: different models have been proposed 
to explain the nature of the X(3872) including S-wave Z}°D*° molecule models [9|, [lOj and 



various diquark-antidiquark models [ill . Il2| . The DD* molecule proposal is motivated by the 



proximity of the X(3872) to the L>°L>*° threshold: m^o+m^.o = 3871.81±0.25 MeV [m,\l£ 
In the molecular model, the X(3872) is a = 1^ state. Some authors have argued [K 
that this model, together with factorization, heavy-quark and isospin symmetries, implies 
that the ratio of B^ X(3872)K0 to B+ X{3872)K+ decays is smaller than 0.1; this 
claim has recently been challenged |l5j. This ratio is expected to be unity for charmonium 
as well as for hybrids (ccg) and glueballs (gg). 

The diquark anti-diquark model of Maiani et al. [ll| predicts that the observed X(3872) 
is one component of a doublet of states. In this model, the X(3872) produced in charged 
B meson decays will have a mass that is different from its counterpart in neutral B meson 
decays by 6M = (7 ± 2)/ cos(20) MeV, where ^ is a mixing angle that is near ±20°. 

In order to test the predictions of these models, we compare branching fraction and 
X(3872) mass measurements in charged and neutral B decays. A previous study per- 



formed by BaBar 16| using 413 fb ^ was not conclusive on these points; this analysis uses a 
larger sample, 605 fb~^ (657 xlO^BB pairs), collected with the Belle detector at the KEKB 



asymmetric-energy e~^e (3.5 GeV on 8 GeV) collider [17| operating at the T(4S') resonance. 



THE BELLE DETECTOR 

The Belle detector is a large-solid-angle magnetic spectrometer that consists of a silicon 
vertex detector (SVD), a 50-layer central drift chamber (CDC), an array of aerogel threshold 
Cherenkov counters (ACC), a barrel-like arrangement of time-of-flight scintillation counters 
(TOF), and an electromagnetic calorimeter (ECL) comprised of CsI(Tl) crystals located 
inside a superconducting solenoid coil that provides a 1.5 T magnetic field. An iron fiux- 
return located outside the coil is instrumented to detect mesons and to identify muons 



(KLM). The detector is described in detail elsewhere [18|. Two different inner detector 
configurations were used. For the first sample of 152 xlO^ BB pairs, a 2.0 cm radius 
beampipe and a 3-layer silicon vertex detector (SVD-I) were used; for the remaining 505 x 10^ 
BB pairs, a 1.5 cm radius beampipe, a 4-layer silicon detector (SVD-II), and a small-cell 
inner drift chamber were used [19|. 
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SELECTION 

Charged tracks are required to originate from the interaction point. A hkehhood ratio 
T^k/tt = I^kI{,I^-k^ I^k)-, where {I^k) is the hkehhood value for the pion (kaon) hypothesis, 
is buih using ACC, TOF and CDC [dE/dx) measurements. For charged kaons, we impose 
T^k/tt > 0.6 that have an 88% efficiency and a 10% efficiency for pions. Kg candidates are 
selected within the 7r"'"7r~ mass range [0.4840, 0.5127] GeV/c^. Requirements on the Kg 
vertex displacement from the interaction point and on the difference between vertex and Kg 
flight directions are applied. This selection is described in detail elsewhere (20| . 

We reconstruct J/ip mesons in the decay channel (/ = e or /i) and include 

bremsstrahlung photons that are within 50 mrad of either the e"*" or e~ tracks (denoted as 
e"'"e~(7)). The invariant mass of the J/ip candidates is required to be within —0.150 GeV/c^ 
< Me+e-(-y) - mj/^ < +0.036 GeV/c^ and -0.060 GeV/c^ < M^+^- - mj/^ < +0.036 
GeV/c^, where mj/^ denotes the J/ijj nominal mass [ij], and Me+e-(7) and M^+^- are 
the reconstructed invariant masses from e^e~(7) and /i"*"//", respectively. The J/ip can- 
didate is then combined with a vr+vr" pair for further analysis: both the X(3872) and 
the iIj{2S), which is used for calibration, decay to this final state. More than one 
J/ipTT^n^ combination may be possible at this stage. An additional cut is applied on 
the M^+^- variable: Mj^+tt- > M {7i~^ n~ J / ip) — {nij/^ + mcut)- For B J/ipn^n^K 
modes {B —> J/ijj7i~^7i~K~^'7i~), we apply the above requirement with mcut = 0.2 GeV 
("^cut = 0.15 GeV). For the former case, this cut corresponds to M^+^- > 389 MeV/c^ for 
the ip{2S) region and M^+^- > 575 MeV/c^ for the X(3872) region and reduces significantly 
the combinatorial background (~ 46% in the charged mode) for a reasonable loss of effi- 
ciency (~ 9%). It also has the property of making the background flat in the M{J /%1)tt^tt~) 
variable. 

To reduce the combinatorial background from e^e" qq continuum events, we require 
i?2 < 0.4 where R2 is the ratio of the second to zeroth normalized Fox- Wolfram moments j^, 
and I cos^^^l < 0.8, where 6b is the polar angle of the B meson momentum in the center-of- 
mass (CM) system, relative to the e"*" beam direction. 

B candidates are obtained by combining a K~^, a Kg or K^tt^ candidate with the 
J/ipTT^TT^ candidate. We select B candidates using two variables: the energy differ- 



ence AE = Eb — -Ebeam = Y,i \J c^Pi + c'^^^ — -^bearrn ^"^^ beam coustraiued mass 



all particles from the B candidate {pi and nii are their CM three-momenta and masses 
respectively) and is the B candidate momentum in the CM frame. If more than one 
candidate is obtained at this stage of the analysis, the candidate with AE closest to zero is 
selected. Only B candidates with |Ai^^| < 30 MeV and Mbc < 5.27 GeV/c^ are considered 
for further analysis. 



From this point onwards, we correct the mass measurement using the known J / ip mass, 
redeflning M{J/ipn^7i~) as M{J/ipn^'7i~)—M{J/ip)+mj/^. The selection cuts isolate a very 
pure sample of B ^ ip{2S)K, ip{2S) J/tpn^n^ decays. These events are used to calibrate 
the M{J/ipTT^7[^) resolution and estimate the systematic uncertainty for the X(3872) mass 
difference. Figure [1] shows the M{J/iIj7t~^tt~) distributions for data near 3686 MeV and 




where the summation is over 



RESULTS 
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FIG. 1: The M( J/V'vr+Tr") distribution for the V'(25') (left) and X(3872) (right) region for charged 
(top) and neutral (bottom) B decays. The curve is the result of the fit described in the text. 

3872 MeV for the charged and neutral B modes. We perform a fit to the M{J /ipTT^TT~) 
distribution to determine the ip{2S) and X(3872) yields, and the signal shape. This fit is 
performed in the region [mj/^ + 0.54, mj/^ + 0.82] GeV/c^. We use the same probability 
density function (PDF) for each signal: a sum of two Gaussians with a common mean. We 
first perform the fit for the charged mode, with the ip{2S) and X(3872) masses and the 
two Gaussian widths as free parameters. The width parameters (mostly determined by the 
large ipi'^S) peak) are then fixed, and we perform the fit for the neutral mode, with only 
the masses as free parameters. This procedure allows a clean comparison of the masses in 
charged and neutral B decay, for both the X(3872) and the ipi'iS) control sample. The 
signal yields are 131.7 ± 15.0 and 27.6 ±6.6 for the X(3872)K+ and the X(3872)ir^ modes 
respectively. The significance is determined from —2 ln(£o/'Cmax) where Cq and £max denote 
the likelihoods returned by the fits with the signal yield fixed at zero and at the fitted value, 
respectively. This quantity should be distributed as x^in^of = 2), as two parameters are 
free for the signal. The calculated significance is then 12. 80" and 5.9(T, respectively. 

Using a Monte Carlo (MC) determined acceptance (e), the results are summarized in 
Table [T] and the ratio of branching fractions can then be calculated: 



TABLE I: X(3872) results obtained in the fit described in the text. 



Mode 


Yield 


<%) 


Significance (a) 


B X ^(X(3872) ^ J/iljTr+ir-) 


B+ ^ 


X{3872)K+ 131.7 ± 15.0 


20.9 


12.8 


(8.10 ±0.92 ±0.66) X lO"*^ 




X{3872)K^ 27.6 ± 6.6 


15.2 


5.9 


(6.65 ±1.63 ±0.55) x lO"*^ 



where we assume the B*^ — >■ X{3872)K^ transition rate to be equal to twice the — »• 
X(3872)i^'^ rate. In this ratio, most of the systematic uncertainties cancel. Therefore only 
the uncertainties due to the T(4S') branching fractions [ij] (2.4%), Monte Carlo statistics 
and MC/data differences are included (Table HTl). The latter source dominates: the differ- 
ences are due to kaon identification (2.2%) and Kg reconstruction efficiency (4.5%). 



TABLE II: Summary of the systematic errors in %. 



Source 


A:(3872)i^+ X{3872)K^ Ratio 




1.4 


1.4 




Secondary B 


1.4 


1.4 


2.4 


MC statistics 


0.2 


0.2 


0.2 


MC decay model 


2.0 


2.0 




Kaon ID 


2.2 




2.2 


Lepton ID 


4.2 


4.2 




Tracking 


6.0 


4.8 


1.2 


Kg reconstruction 




4.5 


4.5 


Total (quadrature) 


8.1 


8.3 


5.7 



As a check, the ratio was estimated for modes with il){2S). The signal yields are 
2916 ± 61 events for B+ ^ ^"(2^)^^+ and 559 ± 25 events for 5° ^ i^{2S)K°s which 
gives Q(^B+^1i!^syK+) ~ ^■'^'^ -'- ^'^^ (stat); the systematic error is the same as that for the 
X(3872) case, ±0.05. The ratio was also estimated using 'ip{2S) — > decays in the same 
dataset, finding q(^b+^1p'(2s]k+) ~ ^-^^ -'- ^-^^ (stat). These two results are in reasonable 



agreement with the ratio calculated from the PDG branching fractions [iJ], 0.96 ± 0.11. 

The difference between the masses in the charged and neutral B modes for X(3872)ii' is 
found to be 6M = Mxk+ - Mxro = (±0.18 ± 0.89) MeV/c^ in data. The same calculation 
is performed for the iIj{2S): from measured masses of (3685.12 ± 0.06) MeV/c^ (charged 
mode) and (3685.23 ± 0.14) MeV/c^ (neutral mode), we find 6M = -(0.11 ±0.15) MeV/c^. 
There is thus no significant evidence of 6M bias; we assign a conservative systematic error 
of ±0.26 MeV/c^ by adding the central value and one-sigma error for the ip{2S), and taking 
the result as a symmetric error. The mass difference between the X(3872) states produced 
in B^ and B^ decay is then 

5M = (±0.18 ± 0.89 ± 0.26) MeV/c^ 
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which is consistent with zero. 

Combining the charged and neutral B samples, we perform a fit to the J/%pTT^Ti~ invari- 
ant mass. We correct the fitted X(3872) mass by the difference between the ^^{28) world 
average ij] and the mass we measure. The corrected mass is: 

M(X(3872)) = (3871.46 ± 0.37 ± 0.07) MeV/c^ 

where the first error is statistical and the second systematic (from the m^(^2S) fit and the 
nominal mass m^(^2S))- 

A similar fit is performed to M{J /ip-K^-K^) for the —>■ {J/tpTT^TT')K^7c~ mode (Fig. [2]). 
The signal yield is 90 ± 19 events. A fit to the Mk-k distribution is then performed to disen- 




FIG. 2: The M{J/ipTr+Tr-) distribution for the -0(25") (left) and X(3872) (right) region for 5° 
J/ipTT^TT^K^TT^ . The curve is the result of the fit described in the text. 

tangle the B° X(3872)i^:*(892)° and 3-body B° X(3872)i^+7r- contributions to the 
final state. We select the events within ±7 MeV around m^(2S) and mx(3872) and fit their 
Mk-k distributions (Fig. [3]). For the B*^ —>■ iI){2S)K^'k^ mode, the PDF for the background 
is the sum of a function {M{Ktx) — uik — m^)° (m^ — m^{2S) ~ Mi^Kn))^, representing phase 
space, and a Breit-Wigner function, to describe the K*{892)^ contribution, obtained from 
M{J/ip7t+n') sidebands (|M( J/V'Tr+Tr") - m^^2S) ± 0.030| < 0.014 GeV). The PDFs for 
the signal are a Breit-Wigner PDF with a free mean and width to represent the i^*(892)° 
component and a histogram PDF obtained from MC to represent the i^2(1430)'' compo- 
nent. The signal yield obtained for the ijj {2 S)K* (892)^ component is 963 ± 44 events and 
corresponds to B {^j {2 S)K* (892)^) = (5.4 ± 0.3 (stat)) x 10"^. This resuh is in reasonable 
agreement with the PDG branching fraction [3], (7.2 ± 0.8) x 10~^. 

For the B^ X{3872)K^7t^ mode, the PDF for the background is the sum of a 
phase space function and a Breit-Wigner PDF obtained from M{J /ipn^n') sidebands 
(|M( J/t/^tt+tt") — mx(3872) ±0.030| < 0.014 GeV). The background yield is fixed from these 
sidebands. The signal is represented by two components: a ii"*(892)° Breit-Wigner PDF and 
a phase space function obtained from MC. The signal yields are 8 ± 10 and 81 ± 20 events, 
respectively. The K*{892)^ contribution is not significant and we set the 90% C.L. limit, 
B{B^ X(3872)ir*(892)°) x B{X{m2) J/^ti+ti-) < 3.4 x 10"^ A product branching 
fraction of B{B'^ X{3872){K+7t-)nr) x i3(X(3872) ^ J/^n+n-) = (8.1±2.0l{i) x 10"*^ 
is also obtained. For the systematic error contributions, in addition to those that enter for 
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FIG. 3: The Kir mass spectrum for the B ip{2S)K7r{Mt) and B X{3872)Ktt (right) can- 
didates. Left: B i;{2S)K*{892f is shown by the dotted red curve, B ij{2S)K^{U30f 
by the dash-dot magenta curve, and the background by the dashed blue curve. Right: B — > 
X (3872) {K+ IT- )nb. is shown by the dash-dot red curve, B X(3872)i^*(892)0 by the dotted 
magenta curve, and the background by the dashed blue curve. 



— > X{3872)K^ (Table [IT|) . we have one more track in the final state, limited statistics 
to fix the background in Mktt (±10%), and possible peaking background contributions in 
M{J/tp7i~^n~), based on a study of the ■ip{2S)K~^7c~ calibration mode, (lio°6%)- 

The result for the X(3872) case is in marked contrast to the ip{2S) case, where the 
non resonant B — > iIj{2S)Ktt component is small and the B^ — > ip [2 S)K* [892)^ and 
5+ —>■ iIj[2S)K~^ branching fractions are of comparable size. K* dominance is also found 
for B J/i/jKtt and xaKu 

In summary, we report the first statistically significant observation of B'^ X [3872) Kg 
decay and measure the ratio of branching fractions to be ^^^JV ^x^^872)k+) ~ 0.82±0.22±0.05, 
consistent with unity. The mass difference between the X states produced in these two decay 
modes is found to be SM = Mxk+-Mxko = (+0.18 ±0.89 ±0.26) MeV/c^, consistent with 
zero. In addition, we search for the X(3872) in the decay 5° X(3872)i\:+7r-, X(3872) 
J/iPn+n-. We measure ^3(5° ^ X{3872){K+7r~)NR) x B{X{3872) J/i^n+n') = (8.1 ± 
2.0lii) X 10-*' and we set the 90% C.L. limit, B{B^ X(3872)ir*(892)°) x i3(X(3872) ^ 
J/^vr+TT-) < 3.4 X 10-6. 
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